Evaluación de fórmulas de transporte de fondo 


en un río de gravas acorazado 


e Raúl López* e Damiáa Vericat * Ramon]. Batalla e 
Universidad de Lleida, Lleida, España 


*Autor de correspondencia 


Resumen 


López, R., Vericat, D., € Batalla, R. J. (marzo-abril, 2015). 
Evaluación de fórmulas de transporte de fondo en un río de 
gravas acorazado. Tecnología y Ciencias del Agua, 6(2), 5-20. 


El curso bajo del río Ebro (NE península Ibérica) presenta un 
cauce con lecho de gravas que discurre aguas abajo de un 
complejo de embalses. En dicho tramo, el lecho experimenta 
ciclos de rotura y restablecimiento de la capa superficial 
acorazada como consecuencia tanto de avenidas naturales 
como de crecidas de mantenimiento. En el citado contexto, 
el objetivo de la presente investigación fue doble. En primer 
lugar, la evaluación de la capacidad predictiva del transporte 
de fondo de tres fórmulas seleccionadas en condiciones de 
coraza rota o alterada. En segundo lugar, el análisis de la 
capacidad de dichas fórmulas para predecir el umbral de 
rotura O alteración de la coraza. Todo ello se fundamentó 
en la comparación entre predicciones y mediciones de carga 
de fondo en el tramo de estudio que abarcaron dos años 
hidrológicos. Conforme a los condicionantes del tramo de 
estudio, las fórmulas finalmente seleccionadas fueron las 
ecuaciones desarrolladas por Parker, Klingeman y McLean 
(1982) (P-K-M), Bathurst (2007) (B) y Recking (2010) (R). De 
acuerdo con los resultados, se recomiendan las fórmulas 
P-K-M y R para la predicción de la carga de fondo cuando 
se rompa o altere la coraza, y se desaconseja la fórmula B 
por su marcada tendencia a la infrapredicción. En cambio, 
para la predicción del umbral de rotura o alteración de la 
coraza se recomiendan las fórmulas B y R. Sin embargo, 
cabe advertir que en este caso ambas fórmulas tendieron a 
predecir umbrales bastante inferiores a los medidos. 


Palabras clave: fórmulas de transporte de fondo, transporte 
de fondo, carga de fondo, transporte de sedimentos, río de 
gravas, lecho acorazado, río Ebro. 


Introducción y objetivos 


La predicción del transporte de fondo adquiere 
particular relevancia en disciplinas diversas: 
ingeniería fluvial, geomorfología fluvial, ecolo- 
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Assessment of Bed Load Transport Formula for an Armoured 
Gravel-Bed River. Water Technology and Sciences (in Spanish), 
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The lower Ebro River (NE Iberian peninsula) has a gravel bed which 
flows downstream from a reservoir complex. In this section of the 
river, the armoured surface of the bed undergoes cycles of breaking 
and reestablishment as a result of natural and control floods. The 
objective of the present work is two-fold. First, to evaluate the capacity 
of three formulas to predict bed load transport under conditions 
in which the armour is broken or disturbed. Second, to analyze 
the ability of the formulas to predict the breakage or disturbance 
threshold of the armour. This was all based on a comparison between 
predictions and measurements of the bed load transport in the study 
section for two hydrological years. Based on the characteristics of 
the study section, the formulas which were finally selected were 
the equations developed by Parker, Klingeman and McLean (1982) 
(P-K-M), Bathurst (2007) (B) and Recking (2010) (R). Given the 
results, the P-K-M and R formulas are recommended to predict 
bed load transport when there is breakage or disturbance of the 
amour, and not formula B because of its considerable tendency to 
underestimate. Formulas B and R are recommended to predict the 
breakage or disruption threshold of the armour. Nevertheless, it is 
worth cautioning that in this case both formulas tended to predict 
thresholds much lower than those measured. 
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gía fluvial, gestión de embalses, ingeniería 
medioambiental, dispersión de contaminantes y 
predicción de desastres naturales (p. ej., Martín- 
Vide, 2013; Re, Kazimierski, € Menéndez, 2014). 
Las presas comportan una discontinuidad en la 
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transferencia de sedimento a través de la red de 
drenaje y modifican el régimen hidrológico de 
la corriente. Si aguas abajo de la presa la ten- 
sión de corte del flujo no supera el valor crítico 
de inicio del movimiento de las partículas de 
mayor tamaño del lecho, pero sí es suficiente 
para transportar los granos más finos, en la 
capa superficial el sedimento deviene más 
grueso y, en consecuencia, puede comenzar la 
formación de una coraza. El grado de acora- 
zamiento del lecho condiciona la distribución 
granulométrica y la magnitud de la carga de 
fondo, por ejemplo, mediante la reducción 
del tamaño de partícula y de la cantidad del 
sedimento transportado (Parker éz Sutherland, 
1990). En general, la coraza resulta estable du- 
rante avenidas de magnitud inferior al caudal 
formativo, pero puede romperse o alterarse al 
paso de crecidas de mayor caudal. Sin embargo, 
la coraza también puede persistir incluso para 
caudales altos, siendo ello concomitante con el 
transporte de partículas de todos los tamaños 
presentes en la capa superficial acorazada 
(Wilcock €: DeTemple, 2005). En ocasiones, la 
coraza se estabiliza (p. ej., condición crítica de 
acorazamiento; Chin, Melville, € Raudkivi, 
1994) y deviene estática, impidiendo que la ero- 
sión progrese. A efectos del presente artículo, el 
término coraza alude a la capa superficial del 
lecho, de carácter semipermanente, formada en 
el curso bajo del río Ebro, tramo que discurre 
aguas abajo de un complejo de embalses. 

En ríos con lecho acorazado, algunos au- 
tores han coincidido en distinguir diferentes 
fases del transporte de fondo en función de la 
movilidad o grado de alteración de la coraza 
(p. ej., Jackson $ Beschta, 1982; Bathurst, 2007; 
Recking, 2010). En caso de que el umbral de 
rotura O alteración de la coraza no haya sido 
alcanzando, el transporte de fondo se califica 
como de fase 1, en la cual la carga de fondo se 
encuentra integrada por sedimento relativa- 
mente fino (proveniente de tramos aguas arri- 
ba o movilizado desde parches de sedimento 
fino presentes en el propio tramo) circulando 
sobre la coraza inmóvil e inalterada. Cuando el 
flujo provoca la rotura o alteración parcial de 


la coraza, se define la fase 2. En la fase 2 del 
transporte pueden participar partículas más 
gruesas constitutivas de la coraza y se inicia la 
incorporación de sedimento fino proveniente 
del material subyacente o subsuperficial a una 
tasa, que es función del grado de alteración 
de la coraza. Por último, una fase 3 puede 
ser concebida una vez que la capacidad del 
flujo logra movilizar todos los tamaños de 
sedimento presentes en el lecho (tanto de la 
capa superficial acorazada como del material 
subyacente). Es preciso advertir que no todos 
los investigadores distinguen entre las fases 2 y 
3, sino que consideran un modelo más simple 
de dos fases. En tal modelo, la fase 1 se caracte- 
riza por un flujo que transporta sedimento fino 
circulando sobre una coraza sin romper. Una 
vez que la coraza se ha roto, se inicia la fase 
2, caracterizada por el transporte de partículas 
gruesas provenientes de la coraza, acompaña- 
das del material fino incorporado desde la capa 
subyacente. En este sentido, la fase 3 ocurre con 
mucha menor frecuencia de modo estricto (es 
decir, movilización de todos los tamaños de 
sedimento presentes en lecho), en especial en 
ríos de montaña de alta pendiente y sedimento 
muy pobremente clasificado (con presencia 
de elementos extremadamente gruesos, con 
frecuencia de origen coluvial), en los que muy 
difícilmente el flujo es competente para trans- 
portar el sedimento de mayor tamaño. 

La mayoría de las fórmulas de carga de 
fondo que de manera tradicional han sido 
usadas fueron parcialmente (y en no pocos 
casos totalmente) desarrolladas a partir de ex- 
perimentos en canales de laboratorio en los que 
se utilizó sedimento de granulometría cuasi- 
uniforme y se daban condiciones de transporte 
en equilibrio (en las que la disponibilidad de 
sedimento no limitaba el caudal sólido) (p. ej., 
Meyer-Peter €: Miiller, 1948; Schoklitsch, 1950). 
Por lo tanto, en su desarrollo o formulación no 
se distinguió entre las diferentes capas o ám- 
bitos granulométricos del material (superficial, 
subsuperficial o la combinación de ambas) 
que pueden ser relevantes en ríos de grava. 
En consecuencia, su aplicación a ríos de grava 
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con lecho acorazado plantea diversas incerti- 
dumbres (en cuanto al diámetro característico 
del sedimento más adecuado a sustituir en las 
fórmulas, así como en lo concerniente al um- 
bral de rotura o alteración de la coraza, y a las 
diversas fases de transporte) y su rendimiento 
es potencialmente limitado. Por consiguiente, 
ha resultado necesario el desarrollo de fórmu- 
las que tuviesen en cuenta de forma expresa el 
efecto del acorazamiento en la carga de fondo 
(p. ej., Parker, 1990; Bathurst, 2007; Recking, 
2010), así como la evaluación o análisis de su 
capacidad predictiva a la hora de aplicarla a un 
río en concreto. 

Los objetivos de esta investigación se 
inscribieron al curso bajo del río Ebro, en un 
tramo ubicado aguas abajo de un complejo de 
embalses, y cuyo lecho experimenta un proceso 
cíclico de rotura y restablecimiento de la cora- 
za. En este sentido, un primer objetivo fue la 
evaluación de la capacidad para la predicción 
del transporte de fondo de una selección de fór- 
mulas en condiciones de coraza rota o alterada 
(es decir, en las fases del transporte 2 o 3). Un 
segundo objetivo consistió en el análisis de la 
capacidad de las ecuaciones seleccionadas para 
predecir el umbral de rotura o alteración de la 
coraza. Dichos análisis se fundamentaron en el 
contraste entre las predicciones de las fórmulas 
y las mediciones de carga de fondo en el citado 
tramo de río. 


Material, métodos y tratamiento de la 
base de datos 


El tramo de estudio y las mediciones en 
campo 


La escorrentía anual de la cuenca del río Ebro 
está muy condicionada por las regiones de 
montaña: el área de montaña representa sola- 
mente alrededor de 30% de la superficie total 
de la cuenca del río Ebro, pero genera cerca 
de 60% de la escorrentía media anual (López 
éz Justribó, 2010). La cuenca del río Ebro drena 
85 530 km? del noreste de la península Ibérica, 
desembocando en el mar Mediterráneo. La 


precipitación media anual de la cuenca es de 
entre 600 y 700 mm. El caudal medio registrado 
en la estación de aforo de la desembocadura es 
de 450 m':s”, lo que equivale a una aportación 
hídrica anual media de 14 300 hm*. Casi 190 
grandes presas regulan 67% (= 7 700 hm?*) de 
la escorrentía anual de la cuenca. El complejo 
de presas más grande de la cuenca se completó 
en 1969 en el curso bajo del río y está integrado 
por tres presas: Mequinenza, Riba-Roja y Flix. 
El conjunto cuenta con una capacidad de 1 750 
hr? (equivalente al 13% de la aportación hídri- 
ca anual media de la cuenca) (Batalla, Gomez, 
£ Kondolf, 2004) (figura 1). 

La hidráulica del flujo y el transporte de 
sedimento fueron medidos de forma regular y 
continua durante las crecidas en la sección de 
muestreo de Móra d'Ebre (SMME) durante el 
periodo 2002-2004. Esta sección de registro se 
ubica 27 km aguas abajo de la presa de Flix y en 
ella el ancho del cauce es de 160 m (figura 1). A 
lo largo de este tramo de registro el río muestra 
un cauce único y de baja sinuosidad, siendo su 
pendiente longitudinal media de 8.5-10*m-m”. 
En dicho tramo, el cauce del río puede conce- 
birse como un canal de ensayos natural, válido 
para el estudio de los cambios en la superficie 
del lecho en función del transporte de fondo. La 
coraza que se ha desarrollado controla la movi- 
lidad del sedimento y la estabilidad del cauce, 
siendo el acorazamiento una expresión de la 
limitación de la disponibilidad del sedimento 
proveniente de los tramos aguas arriba y de la 
frecuencia de ocurrencia de los flujos compe- 
tentes. En este apartado se presenta un breve 
resumen de la metodología de adquisición de 
datos con la que se caracterizó el sedimento del 
lecho del río y se midió el transporte de fondo 
en la SMME. Una descripción pormenorizada 
de dicha metodología puede consultarse en Ve- 
ricat y Batalla (2006), Vericat, Batalla y Garcia 
(2006), y Vericat, Church y Batalla (2006). 

Los años hidrológicos del periodo de estu- 
dio (2002-2003 y 2003-2004) pueden conside- 
rarse medianos en cuanto a caudal medio anual 
(415 m's* en el primero año y 465 m'-s* en el 
segundo). Varias crecidas acaecieron durante el 
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Figura 1. Localización de la zona de estudio en el curso bajo del río Ebro (NE península Ibérica). 


periodo de estudio, midiéndose en la mayoría 
de ellas el transporte de fondo. Algunas fueron 
crecidas de mantenimiento (flushing flows) 
(Batalla €: Vericat, 2009). El caudal máximo re- 
gistrado durante el periodo de estudio alcanzó 
los 2 500 m*-s1 (aconteció en febrero de 2003), 
correspondiéndole un periodo de retorno de 
ocho años, por lo que dicho episodio se con- 
sideró una avenida importante en el contexto 
de la distribución histórica de caudales punta 
de avenida. 

Para el propósito de este estudio, se usó la 
distribución granulométrica del material del 
lecho (tanto del material superficial como del 
subyacente) obtenida en una barra expuesta 
cercana a la SMME. Dicha barra se encuentra a 
menos de 500 m aguas abajo de la SMME (dis- 
tancia equivalente a tres veces el ancho medio 
del cauce) (figura 1). Es el depósito de grava 
accesible más cercano a la sección de medida y 
fue considerado plenamente representativo de 


la distribución granulométrica del sedimento 
activo del río. Por otra parte, el sedimento 
inactivo, el cual fue distinguido por la cubierta 
vegetal, fue evitado, dado que puede guardar 
poca relación con el régimen actual del río. 
El muestreo del material del lecho se llevó a 
cabo en dos ocasiones: a) en el verano de 2002, 
justo antes del comienzo del año hidrológico 
2002-2003 (material del lecho I, MLD; y b) en 
el verano de 2003, justo antes del comienzo del 
año hidrológico 2003-2004 (material del lecho 
IL, MLID. Por lo tanto, pudo establecerse una 
división cronológica de los datos en función 
del material del lecho (DML). En el cuadro 1 se 
muestra el valor de diversos percentiles granu- 
lométricos de ambas muestras o divisiones 
tanto para el material superficial como para el 
subyacente. La capa superficial más gruesa fue 
caracterizada mediante el método de conteo de 
partículas (Wolman, 1954; Rice €: Church, 1996) 
en el caso de MLI. En el caso de MLII, se empleó 
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Cuadro 1. Percentiles granulométricos del material superficial y subyacente del lecho correspondientes al tramo bajo del río 


Ebro (500 m aguas abajo de la sección de muestreo de la carga de fondo). 


Tamaño de partícula Material del lecho I (MLI)* Material del lecho II (MLID? 
D,(mm) Subyacente Superficial Subyacente Superficial 

1D). 10 34 12 19 

D,, 13 39 15 23 

1D, 19 50 21 33 

D,, 26 Ba 26 38 

Da 52 88 48 70 


* MLL: las muestras se tomaron en 2002, al comienzo del año hidrológico 2002-2003. 
» MLII: las muestras se tomaron en 2003, al comienzo del año hidrológico 2003-2004. 


el método de área por peso (Kellerhals éz Bray, 
1971) debido a que se detectó una no desprecia- 
ble proporción de partículas de tamaño inferior 
a 8 mm. El material subyacente fue muestreado 
siguiendo el método volumétrico, previa elimi- 
nación de la capa superficial (Church, McLean, 
éz Wolcott, 1987). 

El número de muestras de transporte de 
fondo utilizables para el propósito de este es- 
tudio resultó de 172 (123 del año 2002-2003 y 49 
del año 2003-2004) (véase el cuadro 2). Aproxi- 
madamente 96% del intervalo del caudal circu- 
lado fue muestreado a efectos de transporte de 
fondo durante el periodo de estudio. La carga 
de fondo se midió empleando un muestreador 
Helley-Smith con boca de 152 mm operado 
desde una grúa. 


El acorazamiento cíclico y la división y 
tratamiento de la base de datos 


El análisis detallado de las características y 
dinámica del acorazamiento del río Ebro en 


el tramo de estudio se describió en Vericat 
et al. (2006a). El lecho del río, en el tramo de 
estudio, experimenta un proceso cíclico de 
incisión y acorazamiento relacionado con la 
magnitud del flujo. El análisis granulométrico 
de las muestras de carga de fondo, junto con 
observaciones de campo durante el periodo 
de estudio revelaron que: a) con posterioridad 
a la primera crecida de diciembre de 2002 la 
coraza persistía; b) las crecidas acontecidas en 
febrero y marzo de 2003 alteraron o rompieron 
la coraza; y c) la coraza puede considerarse 
ya restablecida durante las avenidas de no- 
viembre y diciembre de 2003 (véase el cuadro 
2). Con arreglo a dichas consideraciones se 
dividieron cronológicamente los datos, en fun- 
ción de las condiciones o estado de la coraza 
(DEC), constituyéndose tres conjuntos (cuadro 
2): 1) coraza sin romper (CSR); 2) coraza rota 
o alterada (CR), y 3) coraza restablecida (CRe). 
Asimismo, y en consecuencia, se adoptó una 
hipótesis de asociación entre cada uno de los 
tres conjuntos y la fase de transporte de fondo. 


Cuadro 2. Clasificación de las mediciones del transporte de fondo en función de la fecha de la crecida, de la división del 
material del lecho (DML), de la división del estado de la coraza (DEC) y de la fase de transporte de fondo. 


Crecidas Nb* Na? DML DEC Fase del transporte 
Diciembre, 2002 40 9 MLI CSR Fase 1 
Febrero-marzo, 2003 83 5 MLI CR Fase2 0 3 
Noviembre-diciembre, 2003 49 17 MLII CRe Fase 1 


“Número total de mediciones del transporte de fondo correspondientes al periodo indicado. 
Número de datos del transporte de fondo resultantes de la agrupación de mediciones por clases de caudal líquido y correspondientes al periodo 


indicado. 
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Es decir, se consideró que en los conjuntos CSR 
y CRe, la fase de carga de fondo representativa 
era la 1, y que las fases 2 o 3 caracterizaban el 
transporte de fondo en el conjunto CR (cuadro 
2). Además, nótese la relación resultante entre 
los conjuntos integrantes de las divisiones DEC 
y DML (cuadro 2). 

En la figura 2a se ha representado la rela- 
ción entre las mediciones de caudal líquido 
(Q) y de transporte de fondo (q, expresado 
en unidades de peso) para los tres conjuntos 
delimitados por la condición de acorazamiento 
(DEC). En el cuadro 2 se especifica el número 
de datos que integran cada uno de los tres 
conjuntos. La figura 2a desvela un elevado 
grado de dispersión de los datos; es decir, dado 
un valor de Q, se registra un amplio intervalo 
de variación en el valor de q, (dicho intervalo 
puede alcanzar una amplitud de incluso varios 
órdenes de magnitud). El fenómeno de la alta 
fluctuación temporal de la carga de fondo ha 
sido ampliamente descrito, incluso en con- 
diciones (medias) de flujo constante, como 
consecuencia de estocasticidad, migración de 
formas de fondo, clasificación granulométrica, 
histéresis o limitación de la disponibilidad de 
sedimento (p. ej., Recking, Liébault, Peteuil, €z 
Jolimet, 2012). En la figura 2a, las curvas de re- 
gresión de los conjuntos CR y CSR discurren en 
extremos opuestos, representando diferentes 
condiciones de suministro del material sólido 
de fondo o fases del transporte. Dado un valor 
de caudal, una mayor carga de fondo se espera 
para condiciones CR que para CSR, pues la 
primera representaría la fase de transporte 1 y 
la segunda una fase 2 0 3. La curva de regresión 
de la condición CRe transcurre entre ambos 
extremos, representando una situación inter- 
media. odas estas observaciones deberían 
considerarse con precaución, habida cuenta del 
bajo valor del coeficiente de determinación (R?) 
de las curvas de regresión, en especial la co- 
rrespondiente al conjunto CR. En este sentido, 
una de las técnicas que pueden utilizarse para 
suavizar el efecto de la fluctuación de la carga 
de fondo y la consiguiente dispersión de datos 
es la agrupación de dichos datos en clases de 
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la variable independiente (en este caso, Q). En 
el presente estudio se agruparon los datos de 
los tres conjuntos en clases cuya amplitud de 
intervalo correspondía a 40 m*s”, asignando 
a cada clase un único valor calculado como el 
promedio aritmético de los datos incluidos en 
la clase. En el cuadro 2 se especifica el número 
de datos que constituyeron cada conjunto se- 
gún la agrupación descrita. En la figura 2b se 
observa que la agrupación por clases conllevó 
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Figura 2. Relación entre las mediciones de caudal líquido 
(Q) y de transporte de fondo (q,) para las distintas 
condiciones de acorazamiento: coraza rota (CR), coraza sin 
romper (CSR) y coraza restablecida (CRe). Figura 2a: base 
de datos completa; figura 2b: base de datos agrupada según 
clases de caudal líquido (Q). 
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una atenuación importante de la dispersión y, 
en consecuencia, el valor de R? de las curvas de 
regresión de los tres conjuntos se incrementó de 
manera notable (en especial el de CR). Puede 
apreciarse cómo la posición relativa de las tres 
curvas se mantiene respecto a la de la figura 2a. 


Selección de las fórmulas de transporte de 
fondo en lecho acorazado 


El requisito fundamental para la elección de 
las fórmulas evaluadas fue que hubiesen sido 
desarrolladas ex profeso para la predicción del 
transporte de fondo en ríos de grava, teniendo 
en cuenta de forma explícita el efecto provoca- 
do por el acorazamiento del lecho. Dichas fór- 
mulas debían, al menos, permitir el cálculo de 
la carga de fondo en la fase 2 o 3 del transporte, 
es decir, cuando el flujo fuera capaz de alterar 
o romper la coraza. Asimismo, se impuso como 
requisito que también incorporasen un método 
de cálculo del umbral de rotura o alteración de 
la coraza (es decir, de inicio de la fase 2 o 3). 
Sin embargo, no se consideró preciso que per- 


mitiesen la predicción del transporte de fondo 
en la fase 1, por no ser necesario para cumplir 
con los objetivos de la presente investigación. 
En consecuencia, las fórmulas seleccionadas 
se aplicaron sólo para la condición de coraza 
rota (CR), involucrando tanto al subconjunto 
completo (N = 83, véase el cuadro 2) como al 
subconjunto agrupado (N = 15, véanse los cua- 
dros 2 y 3). Además, como requisito adicional, 
se obviaron de manera expresa aquellas fórmu- 
las basadas en el cálculo de la carga de fondo 
por fracciones granulométricas (p. ej., Parker, 
1990). En la práctica, la aplicación de este tipo 
de ecuaciones es difícil debido, por ejemplo, a 
la frecuente escasez de datos granulométricos 
completos, al mayor esfuerzo de cálculo o al 
mayor número de variables que requieren ser 
medidas en campo (Recking, 2010). 
Dadas expuestas, 
fórmulas finalmente seleccionadas para su 


las restricciones las 
evaluación fueron las desarrolladas por Parker 
et al. (1982), Bathurst (2007) y Recking (2010) 
(referidas en lo sucesivo como P-K-M, B y R, 


respectivamente). En el cuadro 4 se muestran 


Cuadro 3. Valor de las variables hidráulicas para el conjunto agrupado de datos (según clases de caudal líquido y para la 


condición de coraza rota, CR) (véase el cuadro 2). 


ye a Ea a a 
Al 727 144 3.43 1.09 
A2 759 145 3.50 0.81 
A3 810 146 3.62 137 
Al 861 147 SA 0.51 
A5 888 147 3.80 257 
A6 921 148 3.88 2.04 
A7 962 149 3.97 1.33 
A8 993 149 4.05 0.88 
A9 1 044 150 4.17 1.99 
A10 1092 151 4.28 PATO 
A“ 1121 152 4.35 1.34 
A12 1 208 154 AMOO 0.54 
A13 1278 155 4.71 2.00 
A14 1 454 158 DAS. 10.68 
A15 1 493 159 5:22 6.66 


“La sección de estudio puede considerarse hidráulicamente ancha (el valor del radio hidráulico tiende al del calado medio). 
Q es el caudal líquido de la corriente; T, el ancho superficial del flujo; y, la profundidad del flujo, y q,,, es la carga de fondo unitaria medida. 
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Cuadro 4. Selección de fórmulas de transporte de fondo de aplicación a ríos con lecho de grava acorazado. 


Ámbit Int íri 
Referencia Sigla Variables requeridas* Sa 5 miemvalo RS 
empírico o de aplicación 
Parker et al. (1982) P-K-M S, 87 Ps Py Ys Do, Río D,,, < 28 mm 


12 < D,, (mm) < 146 


Bathurst (2007) B [Da ho Dese [Dg (Dl =600 Río 30 < D,, (mm) < 540 
152<D,/D,, < 11 
.02 8 
Recking (2010) R S, 87 Pr PyR, Diy Doo =7600 | Río y laboratorio 0.02 < S(Jo) < 


0.9 < D,, (mm) < 558 


“Variables requeridas para el cálculo de la carga de fondo en las fases de transporte 2 0 3. 


Número de datos involucrados en el desarrollo de la fórmula. 


Ses la pendiente longitudinal del cauce; g, la constante de aceleración de la gravedad; p y p, la densidad del agua y de las partículas de sedimen- 
to, respectivamente; y, profundidad media del flujo en la sección transversal del cauce; R, el radio hidráulico de la sección transversal del cauce; 
Da, y Day el tamaño de partícula de la capa superficial del lecho para el que 84 y el 50% de la muestra de sedimento es inferior, respectivamente, 
y D;,, es el tamaño de partícula del material subyacente o subsuperficial del lecho para el que 50% de la muestra de sedimento es inferior. 


de forma sintética las variables requeridas y el 
intervalo experimental o de aplicación de las 
ecuaciones (nótese que en los tres casos, la base 
empírica incluye datos de ríos). Puede obser- 
varse que tanto el valor de pendiente media 
del cauce como del tamaño de las partículas de 
sedimento del tramo de estudio cumplen con 
los intervalos experimentales de las fórmulas 
evaluadas. La fórmula R permite determinar 
en cuál de las tres fases teóricas se produce el 
transporte de fondo (es decir, fase 1,2 0 3) y pre- 
decir el correspondiente valor de caudal sólido. 
En cambio, la fórmula B sólo predice la carga 
de fondo en el caso de rotura o alteración de la 
coraza (sin distinguir de manera explícita entre 
las fases 2 y 3). Por último, la fórmula P-K-M 
predice el transporte de fondo solamente en 
caso de alteración de la coraza y suponiendo 
que todos los tamaños de partícula inician de 
manera simultánea su movimiento cuando 
acontece la condición crítica de alteración o 
rotura (concepción que guarda mayor relación 
con la definición de la fase 3 del transporte). 


Evaluación estadística de la capacidad 
predictiva de las fórmulas 


La capacidad predictiva de las fórmulas estu- 
diadas se evaluó mediante la comparación de 
los valores medidos (q,,) y predichos (4,,) de 
carga de fondo unitaria (caudal de sedimento 
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expresado en peso y por unidad de ancho, 
en N:s im”). Diversos índices estadísticos y 
métodos gráficos se emplearon para evaluar 
dicha capacidad. Los índices utilizados se 
fundamentan en la relación o razón de dis- 
crepancia (r), definida como el cociente entre 
el valor predicho y el medido (r = q,,/9,,,). El 
intervalo de dicha razón es (0, + oo). En estu- 
dios de transporte de fondo en ríos de grava, 
r puede expandirse en un amplio intervalo de 
valores: frecuentemente de dos o más órdenes 
de magnitud (p. ej., Duan, Chen, éz Scott, 2006; 
Recking, 2010). Por consiguiente, las compa- 
raciones estadísticas entre valores predichos 
y medidos requieren incluir transformaciones 
logarítmicas e índices que resulten menos sen- 
sibles a los valores extremos. 

A continuación se exponen los índices esta- 
dísticos que se emplearon para evaluar el ajus- 
te entre predicciones y medidas de la carga de 
fondo. En primer lugar, se calculó el porcentaje 
de datos cuyo valor de r (q,,/ Q,,) NO SUperase 
una proporción de 2 (0.5<r<2), 5 (0.2<r<5) y 
10 (0.1 <r< 10). Asimismo, se empleó la media 
aritmética de r (mr): 


mr (1/0) 7, (1) 


siendo r, el ¡-ésimo valor de r y N el número de 
datos. Dicho índice varía en el intervalo (0, +00), 
indicando menor discrepancia cuanto más cer- 
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cano sea su valor a 1. También se usó la media 
aritmética del logaritmo de r (mlr): 


N 
mlr =(1/ NY, _logr, (2) 
denotando r, el ¡-ésimo valor de r y N el núme- 
ro de datos. Dicho índice varía en el intervalo 
(ox, + 00), indicando menor discrepancia cuan- 
to más cercano sea su valor a 0. Asimismo, se 
usó una modificación de la media geométrica 
de r (gr) (Habersack €: Laronne, 2002): 
gr= (rre ON (3) 
donde se tomó el valor inverso de r, si r, < 1, 
para asegurar que gr > 1. El índice así definido 
varía en el intervalo (1, + oo), indicando menor 
discrepancia cuanto más cercano sea su valor a 
1. Además, se utilizó una variación ponderada 
del índice gr (gpr) (Habersack € Laronne, 2002): 
epr=(rpyrpy or p yA (4) 
denotando rp un valor de r ponderado por la 
potencia del caudal sólido de fondo medido 
(rp = rísm) y donde se tomó el valor inverso de 
rp, si rp,< 1, para asegurar que gpr > 1. El índice, 
así definido, varía en el intervalo (1, + oo), indi- 
cando menor discrepancia cuanto más cercano 
sea su valor a 1. 

Además de los índices estadísticos expues- 
tos, también se emplearon representaciones 
gráficas de la divergencia entre valores predi- 
chos y medidos, que favorecieron una interpre- 
tación visual del rendimiento de las fórmulas. 
Se representó en coordenadas logarítmicas el 
valor predicho (q,,), en función del valor me- 
dido (q,,,) para cada uno de los datos de carga 
de fondo (diagrama de dispersión). Asimismo, 
se analizó la distribución de la razón de dis- 
crepancia (r) mediante un diagrama de cajas 
representado en escala logarítmica. 

La clasificación de las fórmulas evaluadas 
en función de su rendimiento puede variar de- 
pendiendo de las propiedades estadísticas del 
índice tomado como criterio clasificador. Inves- 
tigaciones previas (p. ej., Barry, Buffington, $ 


King, 2007) hallaron que dada la frecuencia con 
la que los índices utilizados cometen errores 
de sesgo, no existe un índice perfecto para la 
evaluación estadística del rendimiento de las 
ecuaciones. Por consiguiente, lo más recomen- 
dable es el uso de una adecuada combinación 
de varios índices, práctica que se siguió en este 
trabajo. En cualquier caso, para una correcta 
interpretación de los resultados es necesario te- 
ner presente las principales limitaciones de los 
índices utilizados. Por ejemplo, el índice mr es 
más sensible a valores de r mayores que 1 (p. ej., 
un valor de r=10 pesa mucho más en el cálculo 
de mr que un valor de 0.1, pese a que de hecho 
ambos representan una desviación de un orden 
de magnitud respecto al eje de simetría r =1). 
Por el contrario, en el índice mlr, los errores 
de igual magnitud pesan lo mismo, con inde- 
pendencia de cuál sea su posición relativa con 
respecto al eje de simetría logr =0 (p. ej., r =10 
y r = 0.1). Sin embargo, una de sus mayores 
desventajas es que los valores de logr de igual 
magnitud y signo contrario se compensan mu- 
tuamente y dan lugar a mlr = 0. Por lo tanto, 
es más sensible a desviaciones pequeñas, pero 
sin simetría (p. ej., si r, = 1.5 y r, = 2, entonces 
mlr = 0.24), que a grandes desviaciones que 
guarden simetría (p. ej., si r, = 0.01 y r, = 100, 
entonces mlr = 0). Por último, el índice gpr, por 
su propia definición, es más sensible al error 
cometido en la predicción de los valores más 
altos de la carga de fondo medida (q,, ). 


Resultados y discusión 


En los cuadros 5 y 6 se muestran los valores de 
los índices estadísticos para las tres fórmulas 
evaluadas para la condición de coraza rota (CR), 
involucrando al subconjunto completo (N = 83, 
véase el cuadro 2) y al subconjunto agrupado 
(N = 15, véase el cuadro 2), respectivamente. 
En la figura 3 se representó para cada dato la 
relación entre el transporte de fondo medido y 
predicho (según las tres fórmulas evaluadas) 
para la condición de CR y los subconjuntos 
completo (figura 3a) y agrupado (figura 3b), 
respectivamente. Asimismo, en la figura 4, 
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Figura 3. Representación del transporte de fondo predicho por las fórmulas evaluadas en función del transporte de fondo 
medido. Figura 3a: base de datos completa; figura 3b: base de datos agrupada según clases de caudal líquido (Q). Únicamente 
se han representado los datos para los que se predijo, transporte de fondo en fase 2 0 3. Las líneas paralelas a la línea de ajuste 


perfecto (r = 1) corresponden ar =0.1 y r=10. 


para cada fórmula y subconjunto de datos de la 
condición CR (es decir, subconjunto completo y 
agrupado), se representaron los diagramas de 
cajas correspondientes a la distribución de la 
relación de discrepancia (1). 

En primer lugar se analizó la capacidad 
de las fórmulas para la predicción del um- 
bral de rotura o alteración de la coraza. Para 
ello se atendió al valor de los índices de y ndc 


consignados en los cuadros 5 y 6. El índice dc 
se definió como el porcentaje de datos del sub- 
conjunto CR (condición de coraza rota), para 
los que la fórmula predijo rotura de la coraza 
O fases de transporte 2 o 3. En los cuadros 5 y 6 
puede comprobarse que la fórmula P-K-M es la 
que mejor rendimiento mostró, con un 100% en 
ambos subconjuntos (completo y agrupado). 
La fórmula R también alcanzó porcentajes 
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Fórmula de carga de fondo 


Figura 4. Diagramas de cajas de la distribución de la razón 

de discrepancia entre el valor predicho y medido de carga 
de fondo (r). El distintivo (a) indica que la fórmula se ha 
aplicado al conjunto de datos agrupado por clases de Q. 


elevados, entre 80 y 93%, dependiendo del sub- 
conjunto. Por último, el resultado de la fórmula 
B fue más sensible a la agrupación de datos, ya 
que su porcentaje descendió de 81 a 53%. En 
un sentido complementario, el índice ndc se 
definió como el porcentaje de datos de la suma 
de los subconjuntos CSR y CRe para los que la 
fórmula no predijo rotura de la coraza o predijo 
fase de transporte 1. En este caso, el rendimien- 


to de las fórmulas resultó muy inferior al ana- 
lizado antes. En efecto, las fórmulas B y R, que 
mostraron resultados similares, no superaron 
30% en el mejor de los casos, mientras que la 
fórmula P-K-M apenas superó 1%. Por consi- 
guiente, en términos generales, en el tramo de 
estudio, las fórmulas evaluadas predijeron un 
umbral de rotura o alteración de la coraza que 
se encontró por debajo del observado. Es decir, 
predijeron el inicio de la fase de transporte 2 
o 3 antes de que ésta se produjera realmente. 
En función de los resultados expuestos y a 
efectos prácticos del tramo estudiado, lo más 
recomendable sería emplear de forma conjunta 
las fórmulas B y R (evitando, en cualquier caso, 
el uso de la fórmula P-K-M) para predecir el 
umbral de rotura o alteración de la coraza, y 
el inicio de la fase de transporte de fondo 2, 
aunque teniendo en cuenta que muy probable- 
mente dicho umbral resulte demasiado bajo. 
En segundo lugar se analizó la capacidad 
de las fórmulas evaluadas para predecir la 
magnitud de la carga de fondo de los datos 
que integran el conjunto CR tanto para el sub- 
conjunto completo como para el agrupado. Es 
importante señalar respecto a los estadísticos 
consignados en los cuadros 5 y 6, que su valor 
se obtuvo involucrando sólo a aquellos datos 
para los que las fórmulas predijeron rotura o 
alteración de la coraza, es decir, para el por- 


Cuadro 5. Capacidad predictiva de las fórmulas evaluadas respecto a la condición de coraza rota (CR) y subconjunto completo 


(83 datos). Los valores en negrita y subrayados corresponden a los primeros y segundos mejor clasificados, respectivamente, de 


acuerdo con cada índice estadístico. 


ala de ndc> r(0.5-2)% | r(0.2-5)* | r(0.1-10) mr mlr gr gpr 
(%) (%) (Jo) (Jo) (Jo) (0) (0) (9) (09) 

P-K-M 100 1.2 43 70 87 6.27 0.083 34 5.3 
B 8l 30 16 39 64 0.62 -0.685 6.1 170 
R 93 27 47 77 87 12.45 0.388 3.2 3.2 


“Porcentaje de datos correspondientes a la condición de coraza rota (CR), para los que la fórmula predice transporte en fase 2 0 3. 
»Porcentaje de datos correspondientes a las condiciones de coraza sin romper (CSR) y coraza restablecida (CRe), para los que la fórmula predice 


transporte en fase 1. 


“0.5<r<2, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 2. 
40.2 <r<5, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 5. 


*0.1<r<10, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 10. 


"A efectos de comparación y clasificación de la capacidad predictiva de las fórmulas, los valores de nr en el intervalo (0, 1) fueron considerados 


de acuerdo con su recíproco (1/mr). 
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Cuadro 6. Capacidad predictiva de las fórmulas evaluadas respecto a la condición de coraza rota (CR) y subconjunto agrupado 
(15 datos). Los valores en negrita y subrayados corresponden a los primeros y segundos mejor clasificados, respectivamente, de 


acuerdo con cada índice estadístico. 


TA dc" ndc> r(0.5-2)% | r(0.2-5)* | r(0.1-10)* mr! mlr gr gpr 
(Jo) (%) (%) (Jo) (Jo) (E) (53) (E) (E) 

P-K-M 100 0 47 67 93 0.90 -0.33 3.1 6.6 

B 53 19 15 25 50 0.20 -0.88 7.6 1 842 

R 80 15 75 92 100 1.73 0.12 1.7 2.9 


“Porcentaje de datos correspondientes a la condición de coraza rota (CR), para los que la fórmula predice transporte en fase 2 0 3. 
Porcentaje de datos correspondientes a las condiciones de coraza sin romper (CSR) o coraza restablecida (CRe), para los que la fórmula predice 


transporte en fase 1. 


“0.5<r<2, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 2. 

40.2 <r<5, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 5. 

*0.1<»r<10, porcentaje de datos cuyo cociente entre el transporte de fondo predicho y el medido (r) no supera un factor de 10. 

A efectos de comparación y clasificación de la capacidad predictiva de las fórmulas, los valores de mr en el intervalo (0, 1) fueron considerados 


de acuerdo con su recíproco (1/mr). 


centaje de datos especificado en el estadístico 
dc. Lo anterior también se aplicó a la hora de 
representar los gráficos de las figuras 3 y 4. 
Ello debe tenerse presente a la hora de valorar 
la comparación de la capacidad predictiva de 
las tres fórmulas entre sí, ya que el conjunto 
de datos no resultó ser exactamente el mismo. 
Además, como ya se advirtió en el apartado de 
métodos, la clasificación del rendimiento de las 
fórmulas puede variar en función del estadísti- 
co tomado como referencia. 

En términos generales, la capacidad pre- 
dictiva de las fórmulas evaluadas fue relati- 
vamente baja. El valor medio para las tres 
(involucrando al conjunto CR, sin agrupar) del 
porcentaje de datos con una desviación de la 
carga de fondo predicha respecto a la medida 
inferior a un factor de 2 (0.5<r<2),5(0.2<r<5) 
y 10 (0.1 < r < 10) fue de 35, 62 y 79%, respec- 
tivamente. Aunque tal grado de divergencia 
pueda parecer elevado es del mismo orden y 
frecuentemente menor que el hallado en inves- 
tigaciones previas sobre el rendimiento de fór- 
mulas de carga de fondo en ríos de grava. Por 
ejemplo, a modo de muestra de publicaciones 
precedentes, el valor promedio correspondien- 
te a desviaciones inferiores a un factor de 2, 5 
y 10 fue de 19, 44 y 75%, respectivamente, en 
el estudio de Martin (2003); de 11, 25 y 47%, 
respectivamente, en el trabajo de Martin y Ham 
(2005), y de 13, 27 y 34% en la investigación de 
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Recking (2010). Además, debe tenerse en cuen- 
ta que las mediciones de la carga de fondo, en 
el tramo de estudio, fueron cuasi-instantáneas 
(es decir, de corta duración respecto a la dura- 
ción total del episodio de avenida) y se ha en- 
contrado que cuanto mayor es la duración del 
periodo de medida mayor es el rendimiento de 
las fórmulas evaluadas (Recking et al., 2012). 
En la figura 4 se representaron los diagra- 
mas de cajas de la distribución de la razón de 
discrepancia entre el valor de carga de fondo 
predicho por cada fórmula y el medido (r) 
tanto para el subconjunto completo como 
para el agrupado. Al comparar los diagramas 
correspondientes a los dos subconjuntos (com- 
pleto y agrupado) para una misma fórmula, se 
comprobó cómo la eliminación de la fluctua- 
ción por clases (motivada por la agrupación 
de datos) produjo una acusada reducción de 
la dispersión en los valores de la razón de 
discrepancia (r) para las tres ecuaciones. Ello 
fue acompañado de una disminución del valor 
de la mediana de la distribución para las tres 
fórmulas, debido a que la reducción de la dis- 
persión no fue simétrica sino que fue mucho 
más intensa para los valores de r mayores de 1. 
En términos generales, en los cuadros 5 y 6 
se observó que la fórmula B es la que mostró 
una menor capacidad predictiva tanto para el 
subconjunto completo (cuadro 5) como para 
subconjunto agrupado (cuadro 6). La única 
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excepción se dio para el índice mr y el subcon- 
junto completo (cuadro 5). En dicho caso, la 
fórmula B resultó la mejor clasificada debido al 
sesgo del índice mr respecto a los valores de r 
muy superiores a 1, mucho más frecuentes en 
el caso de las fórmulas P-K-M y R (figura 4). La 
fórmula R fue la que mostró un mayor ajuste 
entre predicción y medición, salvo algunas 
pocas excepciones. Dos de las excepciones se 
dieron para el índice mr (cuadros 5 y 6) y fueron 
motivadas por la mayor sensibilidad de dicho 
estadístico con respecto a valores de r muy 
superiores a 1, como se ha expuesto más arriba. 
La otra excepción, que se dio para el subcon- 
junto completo y conforme al valor del índice 
mlr (cuadro 5), se atribuyó a la baja sensibilidad 
de dicho estadístico cuando las desviaciones se 
disponen de forma simétrica (lo que se pone 
de manifiesto en la figura 3a, al comparar los 
diagramas de dispersión de las fórmulas R y 
P-K-M, la ecuación mejor clasificada en este 
caso). La fórmula P-K-M se clasificó mayorita- 
riamente en segunda posición, si bien su valor 
en los estadísticos se aproximó mucho más a la 
fórmula R que a la B (figuras 3 y 4). Asimismo, 
las fórmulas P-K-M y R predijeron con mucho 
menor error los datos de mayor magnitud de 
carga de fondo medida y, por el contrario, para 
los valores más bajos de transporte de fondo 
registrado, el error de predicción fue en com- 
paración mucho más elevado (figura 3). 

Más allá del grado de ajuste concreto entre 
medición y predicción, es de interés conocer si 
una fórmula tiende a la sobrepredicción o la 
infrapredicción del transporte de fondo. En las 
figuras 3 y 4 se observó la acusada tendencia 
de la fórmula B a predecir valores de la carga 
de fondo inferiores a los medidos. Igualmente 
la fórmula P-K-M mostró tendencia a la infra- 
predicción, aunque en un grado mucho más 
moderado. Por último, la fórmula R tendió a 
sobrestimar las mediciones de carga de fondo, 
aunque con un sesgo también comparativa- 
mente leve. Las tendencias descritas lo son 
en términos generales, es decir, abarcando e 
integrando toda la amplitud del intervalo de 
carga de fondo medida. Sin embargo, resulta 


también de utilidad conocer si la tendencia 
varía en función de la magnitud de otras va- 
riables correlacionables con la carga de fondo, 
por ejemplo, el caudal líquido. Para ello, en la 
figura 5 se representaron, en función del cau- 
dal líquido medido (Q), los datos del conjunto 
de condición CR (subconjunto agrupado), su 
curva de regresión y las curvas de regresión 
correspondientes a las tres fórmulas evaluadas. 
La figura 5 develó claramente: la infrapredic- 
ción sistemática de la fórmula B para todo el 
rango de caudal, que para caudales inferiores a 
aproximadamente 1 100 m?/s, la fórmula R es 
la que más se ajustó y que para caudales supe- 
riores a dicho valor se dio poca diferencia entre 
las predicciones de las fórmulas P-K-M y R. 
Cabe advertir que en el desarrollo de la 
fórmula de Recking (2010) se involucraron 46 
datos pertenecientes al subconjunto de condi- 
ción CRe (cuadro 2), tomados de Vericat et al. 
(2006b). Sin embargo, dado que únicamente 
se utilizaron (junto a más de 3 000 datos) para 
verificar el rendimiento de la fórmula y no para 
su ajuste o calibración, ello no compromete la 
independencia de la base datos del presente 
estudio respecto a la fórmula R. Asimismo, im- 
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Figura 5. Relación entre las mediciones de caudal líquido 
(Q) y de transporte de fondo (q,) para el conjunto coraza 
rota (CR) y su curva de regresión (en trazado discontinuo) 
representadas junto a las curvas de regresión de las 
fórmulas P-K-M, B y R para el mismo conjunto de datos. 
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porta destacar que los resultados de la presente 
investigación pueden ser sensibles al periodo 
de muestreo (dos años) y al hecho de que las 
mediciones de carga de fondo se limitaron a un 
único tramo de río. Un periodo de muestreo 
más prolongado hubiese podido revelar dife- 
rencias en la distribución de la carga de fondo 
medida, pudiendo a su vez originar cambios en 
el rendimiento de las fórmulas evaluadas. No 
obstante, algunas de las ventajas de los años 
concretos que se registraron radican en que és- 
tos representan años hidrológicos medios, que 
permitieron observar de manera significativa, 
y poco habitual en ríos de esta magnitud, di- 
ferentes etapas del ciclo de acorazamiento del 
lecho, variando entre un lecho más estabilizado 
debido a la presencia de una coraza bien de- 
sarrollada y otro lecho más móvil a causa del 
efecto de la rotura o alteración de la coraza. 


Conclusiones 


La investigación se considera de aplicación 
práctica directa al tramo de estudio (p. ej., di- 
seño de crecidas de mantenimiento o inyección 
de gravas), aunque la metodología e incluso 
los propios resultados podrían considerarse 
extrapolables a ríos de grava acorazados de 
características hidráulicas, sedimentarias y 
geomorfológicas similares. 

En términos generales, las fórmulas evalua- 
das predijeron umbrales de rotura o alteración 
de la coraza que se encontraban por debajo de 
los observados; es decir, predijeron el inicio de 
la fase 2 0 3 antes de que ésta tuviese lugar. Sin 
embargo, se detectaron importantes diferen- 
cias entre las tres ecuaciones, lo que conlleva 
implicaciones prácticas. En efecto, dado su me- 
jor resultado relativo, se recomienda aplicar las 
ecuaciones B y R como predictores del umbral 
de rotura de la coraza en el tramo de estudio y 
se desaconseja el uso de la ecuación P-K-M. 

Para la predicción de la carga de fondo en las 
fases 2 0 3 del transporte, y después de calcular 
el umbral de rotura o alteración de la coraza de 
acuerdo con las recomendaciones anteriormen- 
te expuestas, se propone la aplicación de las 
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fórmulas P-K-M y R. Ello, teniendo en cuenta 
que, en términos generales, la primera tendió a 
la infrapredicción y que la segunda mostró un 
sesgo hacia la sobrepredicción. Se desaconseja 
el uso de la fórmula B, dado que infrapredijo de 
manera acusada. En cualquier caso, debe tener- 
se presente que el rendimiento de las fórmulas 
hallado en este estudio se encontró dentro del 
orden de magnitud del rendimiento referido en 
investigaciones previas, e incluso fue mayor. 
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